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Abstract 

The model of a classical spacetime foam is considered, which consists 
of static wormholes embedded in Minkowski spacetime. We examine the 
propagation of particles in such a medium and demonstrate that a single 
thin ray undergoes a specific damping in the density of particles depending 
on the traversed path and the distribution of wormholes. The missing 
particles are scattered around the ray. Wormholes was shown to form DM 
halos around point-like sources and, therefore, the correlation predicted 
between the damping and the amount of DM may be used to verify the 
topological nature of Dark Matter. 
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1 Introduction 



The nature of Dark Matter (DM) represents one of the most important and yet 
unsolved problems of the modern astrophysics. Indeed, while the presence of 
DM has long been known |T] and represents a well established fact (e.g., see Refs 
[21 [5] and references therein), there is no common agreement about what DM is. 
In the simplest picture DM represents some non-baryonic particles (predicted 
numerously by particle physics) which should be sufficiently heavy to be cold at 
the moment of recombination and those give the basis to the standard (cold dark 
matter) CDM models. The latter turn out to be very successful in reproducing 
properties of the Universe at very large scales (where perturbations are still 
on the linear stage of the development) which led to a wide-spread optimistic 
believe that non-baryonic particles provide indeed an adequate content of DM. 

However the success of CDM models at very large scales is accompanied with 
a failure at smaller (of the galaxies size) scales. Indeed, cold particles which 
interact only by gravity should necessary form cusps (pdm ~ 1/0 in centers of 
galaxie^l [3] (see also Ref. [5] where the problem of cusps in CDM is discussed 
in more detail), while observations definitely show the cored (pdm ~ const) [B] 
distribution. The only way to destroy the cusp and get the cored distribution 
is to introduce some self-interaction in DM or to consider warm DM. Both 
possibilities are rejected at large scales by observing AT/T spectrum (e.g., see 
Ref. [3] and references therein). By other words DM displays so non-trivial 
properties (it is warm or self-interacting in galaxies, however it was cold at the 
moment of recombination and it is still cold on larger (than galaxies) scales) 
that it is difficult to find particles capable of reconciling such observations. 

These facts support the constant interest to different alternatives of the 
DM hypothesis which interpret the observed discrepancy between luminous and 
gravitational masses as a violation of the law of gravity. Such violations (or 
modifications of general relativity (GR)) have widely been discussed, e.g., see 
Refs. [3 [8]. However, it turns out to be rather difficult to get a modification 
of GR which is flexible enough to reconcile all the variety of the observed DM 
halos. Moreover, the weak lensing observations of a cluster merge in Ref. [9] 
seem to reject most of modifications of GR in which a non-standard gravity 
force scales with baryonic mass. 

The more viable picture of DM phenomena was suggested in Ref. [10] (see 
also references therein) and developed recently in Refs. [ITJ[12]. It is based on 
the fact that on the very early (quantum) stage the Universe should have a foam- 
like topological structure [13] . There are no convincing theoretical arguments of 
why such a foamed structure should decay upon the quantum stage - relics of the 
quantum stage foam might very well survive the cosmological expansion, thus 
creating a certain distribution of wormholcs in the Friedman space. Moreover, 
the inflationary stage in the past [14] should enormously stretch characteristic 

lr The presence of cusps formed by the development of adiabatic perturbations follows 
straightforwardly from the conservation of the circulation theorem in the hydrodynamics. By 
other words the fact that the distribution of DM should have cusps in galaxies is equivalent 
to the fact that DM should represent cold non-baryonic particles. 
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scales of the relic foam. The foam-like structure, in turn, was shown to be flex- 
ible enough to account for the all the variety of DM phenomena [TO] [11] ; for 
parameters of the foam may arbitrary vary in space to produce the observed 
variety of DM halos in galaxies (e.g., the universal rotation curve for spirals 
constructed in Ref. [TS] for the foamed Universe perfectly fits observations). 
Moreover, the topological origin of DM phenomena means that the DM halos 
surrounding point-like sources appear due to the scattering on topological de- 
fects and if a source radiates, such a halo turns out to be luminous too [IT] 
which seems to be the only way to explain naturally the observed absence of 
DM fraction in intracluster gas clouds [5]. 

While the foam-like structure of the Universe is capable of providing a quite 
good description of DM phenomena, it is necessary to look for some indepen- 
dent tests to verify the topological nature of DM. Effects of the spacetime foam 
attract the more increasing attention (e.g., see Refs. [16]-[2Q] and references 
therein). However most of effects considered [H]-[2Q] assume the foamy struc- 
ture at extremely small scales (which correspond to energies higher than 200 
GeV). DM phenomena, however, suggest that the characteristic scale of the 
spacetime foam L (and respectively of wormholes) should be of the galaxy scale, 
e.g., of the order of a few Kpc. The fact that the fundamental length scale for 
the quantum dynamics of spacetime need not be equal to the Planck length was 
also discussed recently in Ref. [T7] , 

In the present paper we consider the propagation of cosmic rays in the foam- 
like Universe. To this end we consider the model of the spacetime foam [12] . 
which consists of a static gas of wormholes embedded in the Minkowski space. 
However contrary to above mentioned papers (e.g., see Ref. [19] where effects of 
cosmic ray interactions in a small-scale foam have been considered)) we assume 
that the characteristic scale of such a foam is of the order of a galaxy size. We 
demonstrate that the scattering on the topological structure is described by a 
specific term in the Boltzmann equation. We show that in a foamed space a 
single thin ray of particles emitted undergoes a specific damping in the density 
of particles depending on the traversed path and the distribution of wormholes, 
while the missing particles are scattered and form a halo around the ray. Such 
halo however has a very low density and is difficult to observe. It turns out 
that the damping traces rather rigidly the amount of wormholes which in the 
foam-like Universe form DM halos in galaxies. Thus, there should exist a rather 
strong correlation between the damping and the distribution of DM in a galaxy 
which presumably can be used to verify the topological nature of Dark matter. 

2 Boltzmann equation 

In the present section for the sake of simplicity we consider the flat Minkowski 
space, while the generalization to the case of Friedman models is straightforward. 
Basic elements of relativistic kinetic theory can be found in standard textbooks, 
e.g., Ref. [21]. Let f (r,p,t) be the number of particles in the interval of the 
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phase space dT = d 3 rd 3 p. This function obeys the equation 

^+r^+p^ = C[f]+a(r,p,t)-\v\ J /3(T,T')f(T')dr' (1) 

where C[f] stands for collisions between particles, a (r,p,t) stands for the rate 
of emission of particles in the phase volume dT, and (3(T, T') describes the 
scattering on wormholes. For the sake of convenience we also distinguished the 
multiplier \v\ — pjra. Our aim is to find an explicit expression for (3(T,T'). 
We consider first a single wormhole, which represents a couple of conjugated 

spheres S± of the radius a and with a distance d = R + — R- between centers of 

spheres. The interior of the spheres is removed and surfaces are glued together. 
The gluing procedure defines the two type of wormholes passable (traversable) 
and impassable. The impassable wormhole appears when before gluing we turn 
out one of surfaces S±. The impassable wormhole works merely as a couple of 
independent spherical mirrors (absolute mirrors, since they reflect gravitons as 
well). The passable wormhole works like a couple conjugated mirrors, so that 
while an incident particle falls on one mirror the reflected particle comes from 
the conjugated mirror. 

Consider an arbitrary point r on the sphere S-, i.e., r € <SL and therefore 

= — R^j = a 2 . The gluing procedure transforms this point into a 

conjugated point f* G 5+ which has the form f = R + + £ + where £ + relates to 
<f_ by some rotation £° = U^^_. Then for the traversable wormhole we find 

J p(T, T')f(T')dT' = (/ - /;) S K+ - a) + (/ - f'_) 6 (£_ - a) , 

where we used the notations £± = r — R± , f± = f (r± ,p±,t), 

f± = R T + U^i ± , 
P± = UT 1 (p-2(pn ± )n ± ), (2) 

and n± = £,±/a. This defines the scattering matrix j3(T, V) for a single wormhole 
in the form 0(T,T') = + (T,T') + 0-(T,T') where 

0± (T, n = 8 K ± - a) [S(r-r')S(p- p') -5(r±-r')6 (p ± - p')] . (3) 

In the case of a spherical mirror (i.e., of the impassable wormhole) this expres- 
sion reduces to the more simple form 

(3{T, V) =5(\r-R\-a)S(r- r') [5 (p -p')-S ( Pl - p')] , (4) 

where p\ = p — 2 (pn) n. 

Let F (R±,a,U) be the density of wormholes with parameters R-, R+, U 
and a, i.e., 

F(R ± ,a,U) = ^2s - R 1 ^ S - R\ ) 5 (a - a n ) S (U - U n ) . (5) 
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Then the total scattering matrix is described by 

/4 ot (T,r') = J [3±(r,T r )F(R±,a,U)d 3 R + d 3 R-dUda. (6) 

We note that the distribution of wormholes has in general quite irregular 
and random behavior and in practical problems it requires some averaging out 
F (R±,a,U), while for a specific astrophysical object (e.g., a galaxy) it may 
possess sufficiently strong fluctuations SF ~ F, 



3 Topological damping of cosmic rays 

In the present section we consider the first terms in the topological scattering 
matrix © and ©. Those terms define the capture of particles by wormholes 
which leads to a specific damping of cosmic rays. Indeed, let us neglect colli- 
sion^] and the topological scattering in (fTJ) and consider trajectories of particles 
x (t) = x (xq,Po, i), p(t) = p (xq,Pq, t). Then we can take variables (xq,Pq, t) as 
new coordinates (instead of {x,p, £)) and the equation (fT]) transforms to 

f t =a(r(t),p(t),t)-\v(t)\f] 1 (r(t))f+\v(t)\ J [3 2 (r,T>)f(r')dT> , (7) 

where 0i describes the capture of particles, while /?2 describes the remission of 
the same particles by wormholes. Now if we consider the case when the source 
a (t) produces a single thin ray and assume that wormholes have isotopic distri- 
bution around the source, then almost all particles captured by wormholes leave 
the ray and will radiate from another regions of space and will have different 
(from the ray) directions. Then in the first order we can neglect the last term 
in r.h.s. of ([7]) and find the solution in the form 

/ = e- T 7, (8) 

where / obey the standard kinetic equation with topological terms omitted 
(i.e., df jdt = dj/dt + rdf /dr+ pdf /dp = a(t)), while the optical depth t (t) 
describes the damping along the ray 

r(t)= /%i(r(i , ))l«(OI* / = /%i(r(s))d S , (9) 

J t JO 

where £ is the coordinate along the ray. 

For astrophysical implications (when the characteristic width of rays L ^> a) 

we can replace S (£± — a) in ([3]) with ira 2 S (^R± — (which means that the 

absorption of particles occurs at the positions R±, i.e., we neglect the throat 
size a). Then, from ^ we find 

/3i (r) = 7T a l S (Rs-r)=n f a 2 n(r, a)da, (10) 

n,s=± J 

2 For the topological damping the absence of collisions is not essential though, since they 
modify merely the function / in JHJ. 
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where n = n + + rt_, and n s (r, a) = J5[R S -?)F (R±,a, U) d 3 R+d 3 R-dU. 
For the sake of simplicity we consider the case when the distribution of worm- 
holes reduces to F (R±,a, U) — g (a) F (R±, U). Then the value (5\ (r) can be 
expressed via the density of wormholes as 

#l (r) = ira?n (r) , (11) 

where a 2 = J a 2 g (a) da and n (r) = n + (r) + n_ (r) is the total density of 
wormholes n± (r) = J2 n $ ( ^± — M • 



4 Topological bias of a point source 

Consider now the case of a stationary point-like source which radiates particles 
in an isotropic way, i.e., a (r,p,t) — X(s)8 (r — fo), where e — \Jp 2 + m 2 and 
A (e) is the distribution of the rate of emission of particles over the momenta. 
Then if we neglect the external force (p = 0), collisions, and the scattering on 
the wormholes the stationary solution to ([1]) is 

fo (r, p) = mX{£ \ 2 6 (cos 6-cosO')5{ V -y') (12) 
p\r- r \ 

where 8, <p define the direction of the vector (r — fo) and 6*', ip' that of p. 

When the density of wormholes is low enough the topological term can be 
accounted for in the next order which defines the topological bias of the source 
a — » a + Sothaio, where the halo density is given by 

Sa halo (T) = \v\ J /3* ot (r,r')/o(r')dr'. 

Such a halo has the two terms Sahaio (r) = 5ai t halo + 5ot2,haio, where the first 
term describes the damping pd|) and the second term defines the remission of 
particles. The exact form of the halo can be found by the image method as in 
Ref. [T2]. Indeed, if we continue the solution to the whole space (we recall that 



the inner region of wormholes 



r-R r ± 



< a n represents the non-physical region 



of space), the wormholes will produce secondary sources of particles. Thus, 
when we neglect the throat size (a <C R±) and assume the isotropic distribution 
over the matrix U, then upon averaging over U every wormhole will radiate in 
the isotropic way which defines the halo as 

Sa 2 ,haio (r,p) = A (e) B 2 (?) , 

where 

B 2 (r)= ^ a 5 (?-!£.), (13) 

?!,S = ± Iff™ — fo 



G 



which defines an additional distribution of particles in the form / (r,p) — 
fo (r,p) + Sf(r,p) 

5f(r,p) = ^± J2 ] ^ ^(cos^-cos^f^.-*/). 

The above expressions can be re- written via the distribution ([5|), e.g., 

B 2 {r)= — — — zN{r,R,a)d 3 Rda, (14) 
J \R-r \ 

where N (r, R, a) =N+ + N_ and N s = J 5 (jl — R- s ^ 5 (r - i? s ) F(R±,a,U) 

d 3 R + d 3 R_ dU (we point out to the obvious relation n(r, a) — J N (r, R, a) d 3 R 
with the distribution n(r,a) in (|10[) ). 

In this manner we see that both functions the damping of cosmic rays (jlOl) 
and the distribution of secondary sources (the halo density) (jT5J) are determined 
via the same function N (r, R,a), i.e., the distribution of wormholes which has 
an irregular (random) behavior. Together with N (r, R, a) functions /3\ (r) and 
£?2 (r) acquire the random character. However, due to the functional dependence 
on the only random function N (r, R, a) such quantities should exhibit a rather 
strong correlation. 

We point out that the interpretatin of the cosmic rays damping possesses 
an ambiguity. For instance a suppression of the cosmic ray flux could be also 
due to other effects, like multiple scattering in the source itself (e.g. see Ref. 
[22] and references therein). Such effects produce analogous correlation between 
the damping and the halo of the secondary sources. Moreover, the halo of the 
secondary sources (fT4"|) is rather difficult to observe; for the brightness of such a 
halo is very low (the intensities of the secondary sources are strongly suppressed 
by the factor a 2 /R 2 , where a is the effective section of the scatterer and R is the 
distance to the scatterer). However, the key point which allows to disentangle 
this specific topological damping from other effects is that the same distribution 
of wormholes defines the distribution of dark matter which we discuss in the next 
section. 



5 Dark matter halos 

As it was demonstrated in Ref. [12) (see also discussions in Refs. (TQl [TT] ) 
the distribution of wormholes defines the density of Dark Matter halos in 
galaxies as well which is much more easier to observe. Indeed, in the presence of 
the gas of wormholes the modification of the Newton's potential was shown to 
be accounted for by the topological bias of sources, i.e., S (i — ro) — > 5 (r — ro) + 
b (r, ro), where the halo density b (r, ro) is determined via the same distribution 
of wormholes ([5]) by expressions analogous to (fT3)l . e.g., see for details Ref. [i~2] . 
The form of the bias function b(r,ro) however admits the direct measurement 
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by observing rotation curves of galaxies (e.g., see Refs. [5] and for the exact 
form of the bias see Refs. [THl HI]). Indeed, in galaxies the topological bias 
relates the densities of dark and luminous matter as 



Pdm(t) = Jb(r- r')p LM {r')d 3 r', (15) 

which for the Fourier transforms takes the form pDAt(k) = b(k)pLM(k). And 
for a point mass it defines the scale-dependent renormalization of the dynamic 
(or the total) mass within the radius R as 

Mtot (R) /M = 1 + 4tt / b (r) r 2 dr. (16) 
Jo 

In observations the amount of DM is defined by the mass-to-luminosity ratio 
M/L. It is assumed that the luminosity traces the distribution of baryons plm 
which is measured by the observing the surface brightness. E.g., spirals can 
be modeled by an infinitely thin disk with surface mass density distribution 
(surface brightness) plm — ae~ r ' RD 8 (z), where Rd is the disc radius (the 
optical radius is R opt = 3.2i?£>). The total dynamic mass is then defined by the 
rotation curve analysis (or by the dispersion of velocities in ellipticals) [6j. 

Observations show that the mass-to-luminosity ratio M tot (r) / L(r) for the 
sphere of the radius r increases with the distance r from the center of the 
galaxy in all galaxies. However if in HSB (high surface brightness) galaxies this 
ratio exceeds slightly the unity within the optical disk M(R opt )/L > 1 which 
means that there is a small amount of DM, in LSB (low surface brightness) 
galaxies such a ratio may reach M(R opt )/L ~ 10 3 . Such a correlation between 
the surface brightness and the amount of DM in galaxies could give an indirect 
evidence for the topological nature of DM; for in accordance to (fTTj) the amount 
of wormholcs defines the damping of cosmic rays and analogously the amount 
of wormholcs defines the amount of dark matter in galaxies [12] . However, the 
basic mechanism which forms such a feature is different (e.g., see Ref. [23] and 
references therein). Indeed in smaller galaxies supernovae are more efficient in 
removing the gas from the central (stellar forming) region of a galaxy than in 
bigger galaxies and this creates the fact that in smaller objects the disc has a 
smaller baryonic density (a lower surface brightness). 

In the general case the relationship between the distribution of dark matter 
and that of wormholcs is rather complicated, e.g., see Ref. |12j . Nevertheless, 
the renormalization of the intensity a point-like source (fT6"| allows us to find a 
rather simple relation between the bias and the density of wormholes on scales 
R ^> d (where d = R+ — R- ) . We stress that the consideration below has a 

rather illustrative (or qualitative) character, while for actual measurements one 
has to use the exact relations in Ref. [T2"] . 

Indeed, the basic effect of a non-trivial topology is that it cuts some portion 
of the volume of the coordinate space. Therefore, the volume of the physically 
admissible region becomes smaller, while the density of particles emitted be- 
comes higher. From the standard flat space standpoint this effectively looks as 
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if the amplitude of a source renormalizes (|16[) . Consider a ball of the radius r 
around a point-like source. E.g., for an isotropic source the number of particles 
emitted in the unit time in the solid angle dfl = r 2 d<f>d cos 9 remains constant dN 
~ fodQ = const, which gives the standard distribution (JT3J) , i.e., fo ~ I/iwr 2 . 

Let us assume that wormholes have an isotropic distribution around the 
source and for the sake of illustration we shall assume that the distribution 
has also the structure F (R±, a,U) = g (a) F (R±, U). Then in the presence of 
wormholes the physical volume is 

V ph (r) = ^7T (r 3 -r>(r)), 

where f2 (r) = 4-7T J a 3 n{r,a)r 2 drda defines the portion of the coordinate 
volume occupied by wormholes within the the radius r and the density of 
wormholes n (r, o) is defined in (fTT)|) . Therefore, the actual value of the sur- 
face which restricts the ball is S p h (r) = -^V p h (r) and we find for the density 
of particles / ~ I j ' S p h (r) which defines the renormalization of the source (|13[) 
I (r) // = A-kt 2 /S p h (r). Absolutely analogously we can use the Gauss diver- 
gency theorem to estimate the renormalization of the gravity source. Indeed, 
the Gauss theorem states that 

/ nVGdS = 4ir [ MS (r) dV = AttM, 

J Jr<B, 
S(R) 

where G is the true Green function (or the actual Newton's potential). Then for 
isotropic distribution of wormholes it defines the normal projection of the force 
as F n (R) = nVG = AttM / S p h (R). This can be rewritten as in the ordinary flat 
space (in terms of the standard Green function Go = — 1/r (i.e., the standard 
Newton's law) and the coordinate surface S COO r — 47ri? 2 ) F n (R) — M' (R) /R 2 , 
where M' (R) /M = ATrR 2 /S ph (R) which defines the bias funct ion in the form 
IflEJ) or 

i r \ 1 d r 2 

iw ^sjw)- (17) 

We stress again that this function admits the direct measurement in galaxies 
[6lll5j. Now by make use of the above expression for V p h (r) we find the behavior 
of the dynamic mass for a point source as 

-^T-- 1 + T^iW) (18) 

where j(r) = ^ir J a 3 n(r,a) da which can be estimated as j(r) ~ |Sf/?i (r). 
Thus, we see that both quantities the damping (i.e., the optical depth r) and 
the amount of DM (the bias b) are expressed via the same function n (r, a) . 

6 Conclusions 

For a homogeneous density of wormholes n (r, a) = n(a) and 0i (r) = /3i = 
const, the damping is determined merely as r (t) = j3i£ where £ is the coor- 
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dinate along the ray. Thus, the damping defines the characteristic scaled L — 
which has the order L ~ a/7 (where 7 = (a/A) 3 , A 3 ~ 1/n is the volume 
per one wormhole, and o is a characteristic size of throats). The parameter 7 
can be extracted from observations of DM in galaxies, while the scale a rep- 
resents here a free parameter which should be fixed from some additional and 
independent considerations. E.g., for the homogeneous distribution of worm- 
holes the value of a defines the amount of dark energy. Indeed, consider one 
wormhole in the Minkowski space. Then the metric can be taken in the form 
(e.g., see Ref. QU) 

ds 2 = dt 2 - f 2 (r) {dr 2 + r 2 sin 2 ddcf 2 + r 2 di? 2 ), 

where / (r) = 1 + 9 (a — r)(% — 1) and 6(x) is the step function^. Both regions 
r > a and r < a represent portions of the ordinary flat Minkowski space and 
therefore the curvature is = 0. However on the boundary r = a it has 
the singularity which defines the scalar curvature as R = —T = -5 (r — a) 
where T stands for the trace of the stress energy tensor which one has to add 
to the Einstein equations to support such a wormhole. It is clear that such a 
source violates the averaged null energy condition, i.e., T — e + 3p < (e.g., 
for the Friedmann space this results in an acceleration of the scale factor ~ t a 
with a = 3 (g^.p) > 1)j i- e -i represents a form of dark energy. Every wormhole 
gives contribution J Tr 2 dr ~ a to the dark energy, while the DE density is e 
= J an(a,r)da ~ 7/a 2 . Thus, the parameter a can in principle be extracted 
from DE density observations. We note however that one has to be careful 
in using such a parameter in galaxies, since in the general case the value a is 
scale dependent (e.g., for the fractal distribution of wormholes the mean value 
is unstable). 

At the optical radius r opt a galaxy can be considered already as a point-like 
source of gravity and, therefore, for estimates we can use (fT5)l instead of (| 1 5|) 
and (fT71) . In HSB galaxies the amount of dark matter within the optical radius 
r opt is rather small M/L > 1 which gives Md yn ( r opt) /M ^1 + 7 (r op t) with 
7 ( r o P t) *C 1 (i.e., A 3 /a 3 » 1) and we can expect the topological damping to be 
negligible (t < 1). In LSB galaxies the mas-to-luminosity ratio may reach M/L 
~ 10 3 which gives 7 (r pt) ~ 1 and we can expect a considerable damping r ~ 1. 
From the qualitative standpoint this feature agrees with the observed correlation 
between the surface brightness and the amount of dark matter in galaxies which 
can be considered as an indirect evidence for the topological nature of Dark 
Matter. The interpretation of such a feature is ambiguous though (e.g., see Ref. 
[25] and for other effects which lead to suppression of the cosmic ray flux see 
Ref. [22j). However we point out that both quantities 7 and r are functions 

3 We point out that such scale has only statistical meaning, since the actual distribution 
of wormholes cannot be utterly homogeneous, otherwise rays could not reach a sufficiently 
remote observer. In particular, there is evidence for the fractal structure of space (e.g., see 
discussions in Refs. |10II15| ) which means that there always exist geodesies along which light 
propagates almost without the scattering. 

4 One can relace / (r) with a smooth function, this however will not change the subsequent 
estimates. 
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of the same random distribution of wormhole n(a, r) and therefore they should 
exhibit a rather strong correlation which may allow to verify the topological 
nature of DM. 

Presumably, astrophysical objects which may also be used to test the topo- 
logical nature of DM are large scale extragalactic relativistic jets in quasars e.g., 
see Ref. [23] and references therein. The smaller jets which are widely observed 
in active galactic nuclei can also be used in LSB galaxies, where the amount of 
DM is considerable. However the crucial step here is the exact knowledge of 
the launching mechanism which may allow to find the discrepancy between the 
predicted profile of a jet and the actually observed one. 
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